Radiative Seesaw Model with Non-zero 6*13 and Warm Dark Matter Scenario 
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With updated experiment result, we study Ma's radiative seesaw model. We try to use fewer 
assumptions to study the feasibility of this model. Both CDM and WDM scenario are taken into 
consideration. In CDM scenario, the fitting of relic abundance shows that a mass spectrum Mi < 
mo < M2 is preferred if theory remains perturbability. By studying the flavour violation process 
H — ► e7 and relic density abundance, we get very strong constraint for the model. In addition, we 
try to explain muon anomaly a M = (g — 2)/2 within the model, but the contribution induced by 
new particles is too small to account for the discrepancy. In WDM scenario, We show that relic 
abundance can be produced by thermal production with subsequent entropy dilution process. The 
entropy dilution can be attained by A^2 decay in which a mass spectrum Mi < M2 < mo is required. 



I. INTRODUCTION 

In physics and cosmology, dark matter and neutrino 
mass are two outstanding problems for Standard Model. 
Many evidences show that our universe contains a great 
amount of non-luminous matter |lj. The observation 
of neutrino oscillation implies non-zero neutrino masses 
[i[ . These results are most important evidences indicat- 
ing the physics beyond Standard Model. 

Small neutrino mass can be generated in several ways. 
Certainly seesaw mechanism is the most attractive one 
for generating neutrino mass [3]. So far we still do not 
know the nature of dark matter. However, we can realize 
some properties of dark matter from observation, such 
as its gravitational effects and that it must be weakly 
interacting. One essential property is that it must be 
stable. To ensure the stability of dark matter, intro- 
ducing a Z2 symmetry is a simple approach which is 
employed in many models. In this scheme, the lightest 
Z2-odd particle cannot decay into other particles, serv- 
ing as dark matter candidate. The best-known example 
is the imparity in supersymmetric models. Since both 
neutrino mass and dark matter problem are two very im- 
portant issues while exploring beyond Standard Model 
regime, we want to investigate the possibility that they 
might be related. This leads to many models trying to 
solve these problems simultaneously. 

One of these is the radiative seesaw model proposed 
by Ma [4[ which gives a simple way to explain both the 
dark matter problem and smallness of neutrino mass. It 
extends Standard Model with a Z2 symmetry, and uses 
one loop effect to generate small neutrino masses. In ad- 
dition, one scalar doublet and three Majorana fermions 
are introduced. So this model is a modest extension for 
getting dark matter and neutrino mass. 

Some works have already been done along this line. 
Sierra et al |E[ consider a specific mass spectrum, and 
discuss the dark matter phenomena. They deduce some 
collider signatures which can be tested in LHC. They 
also discuss warm dark matter (WDM) scenario, but 
fail to explain the relic abundance. Suematsu et al Q 
discuss cold dark matter (CDM) abundance by study- 
ing the structure of mass mixing matrix and coannihi- 



lation effect. In addition, they take into account the 
\i — >• cy process to get constraint. However, their anal- 
ysis assumes that #13 = which has been found to be 
non-zero in the recent Daya-bay's result. In this paper, 
we want to re-examine Ma's model in light of this new 
experimental result. 

Recently, many experiments on neutrino physics have 
great progress. A large amount of data has been accu- 
mulated so that we are able to determine mixing angle of 
mass matrix and mass difference more precisely. In last 
two years, experiments including T2K [7], Daya-bay @, 
and RENO [9[ present new data about #13. These re- 
sults confirm non-zero 6 is and one CP-violation phase 
needs to be studied. With these information, we can 
now use neutrino mass mixing matrix to study Ma's ra- 
diative seesaw model more closely. 

Recently, WDM scenario gains considerable attention 
Some literatures argue that WDM scenario 
might be able to overcome some problems that typi- 
cal CDM have met in structure formation and galaxy 
halos [l3l-[l5|. As a result, we shall take WDM scenario 
into consideration as well. WDM can be attained by 
several different mechanisms. For example, Dodelson 
and Widrow [16] propose that WDM can be produced 
by active-sterile neutrino oscillation during big bang nu- 
cleosynthesis (BBN) epoch. Shi and Fuller [17J propose 
that WDM can be accomplished by a non-thermal res- 
onant production. In this model, there is no mixing 
with active neutrino. Hence, Dodelson and Widow's 
approach cannot apply to this case. Since the Yukawa 
coupling of the sterile neutrinos in this model is not 
necessarily small, we can realize the WDM scenario 
of Ma's model by thermal production [l8|. Some re- 
searches show that correct dark matter relic density can 
be obtained by thermal overproduction with a subse- 
quent entropy dilution from some decay channel jll| . 

In Ref. [l9| , the WDM of radiative seesaw model from 
thermal production has been studied. However, its en- 
tropy dilution is accomplished by adding another real 
scalar singlet which couples to Higgs doublet. Actually 
extra particle for providing reheating process is not re- 
ally needed. We try to realize required entropy dilution 
within this model, without any additional particle. 
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In this paper, we focus on the Ma's radiative seesaw 
model, trying to study its phenomena. With updated 
experimental result, we can reconstruct the mass ma- 
trix. Different from previous works, we use another 
method to analyze mass mixing matrix. Both CDM 
and WDM would be taken into account. By studying 
the flavour violation process, we would put constraint 
on the parameter space of this model. In addition, we 
pay attention to muon anomaly and try to explain the 
experimental result within the model. 

This paper is arranged as follows. In next Section, we 
briefly introduce the radiative seesaw model. In Section 
III, we study the relic abundance. In Section IV and V, 
we study the phenomena of flavour violation process and 
muon anomaly. Finally we give a summary in Section 
VI. 



II. MODEL 

Ernest Ma's radiative seesaw model Q extends Stan- 
dard Model symmetry into : SU(2) x U(l) x Z 2 where 
a discrete Z 2 symmetry is introduced. In addition, one 
scalar SU(2) doublet 77 = (?7 + ,?7 ) and three Majorana 
fermion Ni (i = 1,2,3) are added. All the new par- 
ticles are assigned as Z2-odd. Therefore, both neutral 
particle 77 and the lightest Majorana fermion Ni could 
serve as the dark matter candidate. Here we suppose 
the dark matter is Ni. With these additional particles, 
the Yukawa coupling of leptons and scalar sector have 
to be modified. The Yukawa coupling of leptons is given 

by 

£y = fap{v a 4> + + L<P°)h + Ki{v a nf - l a r) + )Ni + h.c. 

(1) 

where a, f3 = e,[i,r, and h a i are unknown parameters. 
Majorana mass term of Ni is given by 

^MiN^Ni + h.c. (2) 

where Mi < M 2 < M 3 . 

The Higgs potential of this model is given by 

V{$,rj) =/i 2 $t$ + ulrfrj + Ai($ f $) 2 + \ 2 (rfv) 2 

+ \ 3 {&$){rfv) + A 4 ($ t 7?)(7? t $) 

+ ^A 5 [($V + h . c .] (3) 

where A5 is taken to be real without losing any general- 
ity. In this potential, we will take fi 2 > so that only 
get a non-zero vacuum expectation value (0°) = v 
which causes spontaneous symmetry breaking (S SB) and 
rf will not develop vacuum expectation value to keep 
the Z 2 symmetry unbroken. The neutral component of 
77 can be separated into real part and imaginary part, i.e. 
V — (Vr + ^ 7 7?)/v / 2- From SSB, the masses of physical 



scalar bosons would be split, and the values are given 
by 

m 2 ^) = fi 2 + \ 3 v 2 

m 2 (rg) =fi 2 + X 3 v 2 + (A 4 + \$)v 2 (4) 
™ 2 0?i°) = lA + W + (A 4 - A 5 )v 2 . 

Since A5 is related to the active neutrino masses, it' s 
assumed to be very small. ( The relationship between 
A5 and neutrino masses will be discussed later. ) Hence, 
for simplicity, we assume the masses of rf R and rjj to be 
degenerate, and their mass is given by m§ = fi 2 + (A3 + 
\ 4 )v 2 . 

s / 
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FIG. 1. One- loop generation of Majorana neutrino mass. 

Due to Z2 symmetry, Dirac mass of neutrinos is for- 
bidden. Neutrino masses can only be generated by one- 
loop radiative effect involving rf and Ni as showed in 
FIG. [H And the generated mass mixing matrix of neu- 
trino is given by [6[ 
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Map =^2h ai hpiAi (5) 

i=l 

where is given by 

X 5 v 2 /MA 

^)~(l + T ^^). (7) 

Here we see that we can take A5 to be small to get 
small neutrinos masses. In this model, there are nine 
Yukawa coupling constants, h a i. In order to make our 
analysis simpler, we will reduce the number of unknown 
parameters. Since recent experiments of neutrino oscil- 
lation measure the mixing angles and mass differences 
more accurately, we can use these results to fix some of 
our parameters. 

In general, this mass mixing matrix in eq. (j5j) can be 
diagonalized by a Pontecorvo-Maki-Nakagawa-Sakata 
Matrix (PMNS) which is of the form 
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C12C13 Si 2 Ci3 <§13e ^ 



U = -S12C23 - C 12 S23Sl3e ld C12C23 - Sl2S23Sl3e td S23C13 (8) 



-S12S23 - C 1 2C23S 13 e tS -C12S23 - S 12 C23Sl3e lS C23C13 



where Qj = cos^j , Sij 
phase. The diagonal mass matrix is given by 



sin^j and 5 is the CP violation 



U T MU = diag(rai, ?7i2, 777,3) 



(9) 



where mi, 777-2, and 777,3 are the masses of neutrinos. 
In general, there could be two more Majorana phases 
in this formula. Nevertheless, in order to simplify our 
analysis, we take them to be zero. 

According to the most updated fit [20|, the mixing 



angles are sin #12 = 3.07 x 10 



sin 2 #23 



3.98 x 10" 



and sin #13 = 2.45 x 10 . The CP violation phase 
S = 0.897T. The mass square differences are Am^ = 
7.54 x 10- 5 (eV) 2 and Am 2 31 = 2.43 x 10- 3 (eV) 2 . Here 
the normal neutrino mass hierarchy is assumed, i.e. 
mi < 777,2 < 777,3. From the cosmological observation 
of WMAP, there is an upper bound for total neutrino 
mass, ^nrty < 0.58 eV (95% CL) [l|. Here we assume 
the extreme situation: 777,1 = 0, 777,2 = 8.68 x 10 -3 eV, 
and 777,3 = 0.049 eV. By using these data and some as- 
sumptions, we can construct the mass mixing matrix 
as 



M 



1.18 + 0.37* 
0.96 + 0.16* 
2.58 + 0.50* 



0.96 

1.24- 

1.40 



-0.16* 
0.005 
-0.14* 




As we know, the value of A^ depends on the ratio, 
(Mi /mo). If M3 is much larger than Mi, M 2 and mo, A3 
would be negligible compared to Ai and A2. It means 
that we can neglect the contribution from A 3 and h a s 
in this stage. In this way, we reduce the number of 
parameters that we have to handle. 
We define a real positive number a, 



A2 

Ai 



and five complex parameters, 

ft/xl 



ai 



CL2 



hrl 



e2 



7/i2 



h. 



K2 



(11) 



(12) 



(13) 



By using the values of mixing matrix, we can get a\ = 
0.70 - 0.72* and a 2 = 2.27 + 0.61*, which are useful 
in dealing with the relic density. 61, 62, and 63 can- 
not be exactly solved, depending on a. Since the mix- 
ings in neutrino mass matrix are large, it's reasonable to 
take coupling constants to be of roughly of same order 
of magnitude. Although we take the extreme situation 
that mi = 0, we can also consider other situations and 
get similar results. 



III. RELIC ABUNDANCE 

A. Cold dark matter scenario 

In order to explain the dark matter problem, it's nec- 
essary for a model to have right amount of relic abun- 
dance. Here, we first discuss the CDM scenario, and we 
will also study the WDM in the next subsection. At 
early universe, the temperature is very high. Various 
particles are produced and annihilate, reaching equilib- 
rium quickly. As universe expands, the temperature 
would drop, and particles would decouple from each 
other. Since Ni is absolutely stable, the abundance from 
thermal production would remain as dark matter relic. 
To discuss the magnitude of density, we need to study 
the Boltzmann equation given by [2l| 



dn 
~dt 



+ (3H)n = -(*\v\}(n 2 -n 2 eq ) 



(14) 



where n is the number density of dark matter, H is the 
Hubble constant, and a is the annihilation cross-section 
of dark matter particles, v indicates the relative veloc- 
ity of two annihilation particles. n eq means the number 
density in equilibrium which yields to Boltzmann distri- 
bution. 

In CDM scenario, dark matter particles are non- 
relativistic at decoupling. In general, the annihilation 
cross-section takes on this form: (a \v\) = a e ff + 6 e ff ( I v\ 2 ) 
where terms of higher order in |*7| are neglected. The 
approximate solution for Boltzmann equation is given 
by HU 



^DM 




M p i(a e ff + 3b e f?/x F )p c 



Mi 
(15) 

where Tp is freeze-out temperature at decoupling, g* is 
the number of degrees of freedom(DOF) at freezing-out, 
S is entropy density, M p \ is Plank mass, and p c is critical 
energy density of the universe,. 

In the nonrelativistic limit, the annihilation cross- 
section of N1N1 —> l a lp multiplied by the relative ve- 
locity is given by [6| 



a v 



1 M?(Mj + m*) 2 
48tt (Ml+mlY 11 



a(3 



h.c. 



(16) 



There are two Feynman diagrams (the lowest order) 
which contribute to this cross-section. Because of the 
Fermi statistics there is a minus sign between these dia- 
grams and as a result the s-wave parts of the scattering 
cancel out. Only p-wave part remains. With the value of 
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M x [GeV] 

FIG. 2. Relationship \h e i \ - Mi where R = mo /Mi. This 
plot is based on the non-baryon relic density Qdm^ 2 = 0.110. 
Different values of R = 2, 5, 10, 20 are considered. 



cross-section, we can get the relationship between \h e ±\ 
and Mi. Using current data from several observations, 
the non-baryon relic density is Qdm^ 2 = 0.110 ± 0.006 
with 68 % C.L. [22]. Fitting the relic abundance, we 
get the fit in FIG. H In FIG. El we plot the relation be- 
tween \h e i\ and Mi for different values of R = m§/M\. 
Cases R = 2, 5, 10, 20 are considered. In some re- 
gion of Mi, the required value of \h e \ \ is larger than 1. 
It means that we have to deal with a non-perturbative 
interaction. This kind of situation is even worse for 
larger R cases. We can see that in R = 10 case, the 
coupling is beyond perturbative regime for Mi > 60 
GeV. In R = 20 case, the coupling is beyond pertur- 
bative regime for Mi > 20 GeV, and our calculation of 
cross-section is probably not reliable. Nevertheless, we 
can conclude that a coupling larger than 1 is needed in 
these regions. If we suppose the Yukawa interaction is 
perturbable, a scenario R < 10 is preferred. Since there 
is no reason to require Mi and M2 to be of the same 
order, M2 > 10Mi is an acceptable choice. Hence, a 
mass spectrum Mi < mo < M2 is preferred. 

In following sections, we will take some phenomena 
into consideration, and constrain the relationship we get 
from relic density fitting. 



B. Warm dark matter scenario 

We now study the WDM scenario, where the dark 
matter particle has a mass of few keV. Since the mass 
is quite small, the particle would still be relativistic at 
freeze-out point. The quantity Y = n/S indicating the 
ratio between number density n and entropy density S 
is given by (2l| 



Y Nl 



45C(3) £ eff 



2tt 4 



9* 



(17) 



where g e R is the effective number of DOF of N\ and 
g* s ~ 10 2 is the the number of DOF for the entropy. To 
get the correct non-baryon abundance, Yjy 1 must satisfy 
the condition Yn 1 M\S = ^dmPc where p c is the critical 



density and ^dm is the density fraction of non-baryon 
dark matter. The required value for g e ^ is given by 



„ ^ AkeV\ 

9eS = 0.144 (—y 



(18) 



For sterile neutrinos, there exists a lower bound around 
5.6 keV from high-redshift Lyman-a forest observations 
[23| . If Mi ~ 10 keV, g e ff is about 2-order larger than 
the required value. Thus, this value is too large for 
correctly accounting for current relic abundance. For- 
tunately, It's possible to use an out-of-equilibrium de- 
cay to get around the problem. In this case, we can 
have other particle which has been decoupled and would 
decay into other relativistic particles, thermalizing the 
universe. Under the decay process, extra entropy would 
be produced. Therefore, during the process, entropy 
density would increase: 



S -> 7S , Y 



Y 

7 



(19) 



where 7 is the factor which labels the entropy change. 
And the value of Y also changes under the decay pro- 
cess, decreasing by a factor. 

If Mi < mo < M2, the major decay channel of N 2 



would be N 2 
are given by 



rfv and N 2 — »• r] + l . The decay widths 



T{N 2 -+ rfv) = T(N 2 -+ rfP) c, Ml 



£ 1 



(20) 



The scalar particles are also unstable, decaying into N±. 
The decay widths are given by 



r(r,°^iV lt ,)=r(r / ± ^iV 1 Z ± )^^ Yl 



n (21) 

With these decay rates, following Ref. [18j, we can es- 
timate the value of 7, 



7 = 1-83 g, 



./4 M 2 Y N2 



Ttm^ 



(22) 



where F is the total decay width of For M2 ^ 1 GeV, 
the Yukawa coupling constant is h a i ~ 10 -12 . However, 
since A^'s have no mixing with active neutrinos, WDM 
may not be attained by thermal production with such 
small couplings [lfj. Moreover, if we have such small 
couplings, we could have neutrino with around 1 eV 
Dirac mass in SM scenario by adding a right-handed 
Dirac neutrino. Instead of A^, rj decay is also a can- 
didate for serving the reheating process. But the de- 
cay of r] has been discussed before [24|, and it has the 
same problem. A very small Yukawa coupling < 10 -9 is 
needed for correct relic density. As a result, it motivates 
us to consider other scenarios. For Mi < M2 < mo, A^ 
decays into N\ and two leptons via an virtual rj. We 
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FIG. 3. Relationship between £ and Mi with strong and 
weak BBN bound. Different values of r = mo/M2 are taken 
into account. We consider the cases of 7 = 100 and 7 = 1000. 



can estimate the decay width of this process 

a, (3 



IV. LEPTON FLAVOUR VIOLATION 

As we mentioned in the introduction, this model 
is motivated by the dark matter and neutrino mass 
problems. The existence of dark matter candidate is 
achieved by adding a Z 2 symmetry and neutrino mass 
is completed through radiative effect contributed by ad- 
ditional particles Ni and 77. However, these extra pat- 
icles not only provide the effective neutrino masses , 
but also contribute to some other effects. One signif- 
icant effect is lepton flavour violation process. From 
this kind of process, we are able to get constraints on 
the parameter space of h a i. Due to the smallness of 
neutrino mass, the decay widths of lepton flavour vio- 
lating processes coming from neutrino masses are much 
smaller than present experiment limit, e.g. \i — >• ej, 
Br(/i —> cy) < 10 -40 . However, the additional parti- 
cles in this new model would contribute significantly to 
lepton flavour violating decay l a — ^ Ipj. Therefore, the 
experiment data would put constraint on the model by 
studying these processes. These processes are showed 
by FIG.[H The branch ratio of l a — > ^7 is given by [27j 



^2 \ h z<xhip\ 

a, (3 



Ml 



11227r 3 m^ 



(23) 



where the mass of 77 ± are assumed to be mo for con- 
venience. This decay channel shows great potential to 
dilute the relic density and give the correct value of Yn x • 
Here, one essential constraint from BBN has to be taken 
into account. In order to preserve the normal BBN, the 
lifetime of N 2 should not be too long. In this way, the 
decay process could end before nucleons freeze out. We 
can estimate the reheating temperature to be [18[, 



rh = - 



45 

27T#* 



1/4 



(24) 



where T is the decay width of processes N 2 —> N\ll. 
T r h must be larger than around 0.7 MeV [25[ (weak 
BBN bound) to 4.0 MeV [26] (strong BBN bound). The 
relationship between coupling constant C = S l^2a^i/3| 2 
and M 2 is plotted in FIG.[3]with the lower bounds. Here 
r indicates the ratio (mo/M 2 ), and different values of r 
are considered in the plots. From the plots, we can see 
that only a small region has been excluded. Except for 
M 2 being around a few GeV, most area of M 2 is safe 
for both weak and strong BBN bounds. It shows that 
considering N± as a WDM candidate is a viable scenario. 




Ni 



FIG. 4. Flavour violation process induced by r\ and Ni. 



Br(7 a -s-Z/37) 



647rG|,m; 



(25) 



where a em = e 2 /47r is the electromagnetic fine structure 
constant, Gf is the Fermi constant , and the function 
F(x) is defined by 



F(x) 



1 - 6x + 3x 2 + 2x 3 - 6x 2 lnx 
6(1 -x) 4 ' 



(26) 



This function decreases as x increases. As we men- 
tioned in CDM abundance fitting, the mass spectrum 
Mi < mo < M 2 is preferred to stay within validity of 
perturbative calculation. So in this mass hierarchy, Ni 
play the most significant role in these processes. As a 
result, it gives constraints for the Yukawa couplings of 
iVi, h a \. According to the latest experiment results, 
there are upper bounds for branch ratios Br(r urA 
< 4.5 x 10~ 8 [iH and Br(/i -+ ej) < 2.4 x 10" 12 Q. 
Between these two bounds, the process /j, ej provides 
stronger constraint for our model. Hence, we use this 
upper bound to cut the curve that we get from CDM 
relic abundance fitting. FIG. [5] shows the result after 
using the upper bounds from decay process \i — >• ej. As 
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FIG. 5. The fitting curve from CDM relic abundance fitting 
is presented by red solid lines. Upper bound of \h e \\ given 
by the constraint Br(/x — > ej) < 2.4 x 10 -12 is presented by 
dash lines. Different values of R = 2, 5, 10 are considered. 



in the section of relic abundance, we consider different 
values of ratio R. For R = 2 case, the region of Mi 
larger than 170 GeV is ruled out. For R = 5 case, only 
mass smaller than 20 GeV is allowed For R = 10 case, 
the coupling constant is larger than 1 in some region, 
i.e. the non-perturbative Yukawa interaction regime. 
Besides, almost all the parameter region is excluded by 
the requirement from \i — >• ej. It implies that in CDM 
scenario mo cannot be too larger than M±. This re- 
sult is consistent with the requirement of the theory 
being perturbative. It also suggests the choice of the 
mass spectrum Mi < mo < M2. Although our result 
is based on the analysis of mixing matrix in Section II, 
there would be no big difference because a\,ai are still 
of order one for such a highly mixed mass matrix. By 
analyzing these plots, we found that it also gives a con- 
straint on mo that it must be not larger than 300 GeV. 



The well measured value of branch ratio provides a very 
strong constraint on this model. Even in the region that 
is not excluded, the fitting curve is quite close to the up- 
per limit. It means that if the measurement of \i — >• 
can be improved, it would be possible to detect the ef- 
fect of flavour violation decay induced by Ni and rj. On 
the other hand, if we still do not see any discrepancy be- 
yond standard model in further experiment with better 
precisions, more region would has to be ruled out and 
endanger the validity of this model. 

V. MUON ANOMALOUS MAGNETIC 
MOMENT 

The additional particles, Ni and 77, also contribute 
to the magnetic moment of leptons. Similar with FIG. 
2J but the flavour of incoming and out-going leptons 
are now the same, i.e. a = f3. Among all kinds of 
leptons, the muon magnetic moment shows some inter- 
esting phenomenon. In 2001, E821 experiment at BNL 
found that there is discrepancy between the experimen- 
tal value and prediction of muon anomaly a M = (g — 2)/2 
from the Standard Model (30[. It's possible to explain 
this deviation by new contribution. Thus, here we pay 
attention to the anomalous magnetic moment of muon. 
The muon anomalous magnetic moment contributed by 
additional particles is given by [23] 

A «M = T7T^£lM 2 ^(^f) (27) 
167r 2 mQ ^ \ m o J 

where m M is the mass of muon. For muon anomaly, the 
discrepancy between experiment and theoretical predic- 
tion based on Standard Model is given by a^ xp — a^ M = 
(26.1 ±8.0) x 10- 10 [p. FIG.[6]shows the fitting curve 
of relationship, | h e \ | to Mi, and the required value to 
account for the discrepancy of muon anomaly. From 
these figures, we can see that large coupling constants 
are needed for explaining muon anomaly. The values are 
very close to 1, i.e. the perturbative limit. For R = 2 
case, the coupling constant is larger than 1 for Mi larger 
than 30 GeV. For R = 5 and R = 10, the coupling is 
below perturbative limit only while if Mi is a few GeV. 
The additional loop correction contributed by Ni and rj 
cannot fully explain the observed muon anomaly. 

VI. SUMMARY 

With the latest experimental data and most updated 
fit, we study Ma's radiative seesaw model with fewer 
assumptions and different scenarios. 

In CDM scenario, we have done a more refined anal- 
ysis than the previous work with the new information 
of #13 and mass mixing matrix . Furthermore, if we re- 
quire the theory to remain in the perturbative regime, 
the mass spectrum Mi < mo < M2 is preferred. By 
studying the flavour violation decay \i — » ej, a very 
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FIG. 6. The fitting curves from relic density is presented by 
red lines. And required values for explaining the observed 
muon anomaly is presented by blue lines. Different values of 
R= 2, 5, 10 are considered. 



strong constraint is obtained. The ratio R = mo /Mi 
is favoured not to be larger than 10, and it show a 
upper bound roughly 300 GeV for mo. Compared to 
Ref. [6], the constraint for Mi is more precise, re- 



stricting it to be in the region below 200 GeV. This 
result is consistent with the preferred mass spectrum 
Mi < mo < M.2- In Ref. [5], it study a specific mass 
spectrum Mi <C M2 < M3 < mo- However, according 
to our result, it seems not viable in CDM scenario. We 
try to use additional contribution from new particles 
to explain observed muon anomaly. Unfortunately, the 
muon anomalous magnetic moment induced by N\ and 
r] is too small to fully account for the discrepancy. 

In WDM scenario, we use a new approach. We show 
that WDM relic can be attained by thermal production 
with a subsequent entropy dilution process. The en- 
tropy dilution can be accomplished by N2 decay within 
the model, and no modification for this model is needed. 
In this scenario, the mass spectrum Mi < M2 < mo is 
preferred for providing a required decay channel N2 —> 
Nil! via, an intermediate virtual rj. Different from previ- 
ous works, extremely small couplings are not required. 

Direct detection of dark matter is also an important 
issue and has potential to give some restrictions for 
this model. But, in this paper, we do not discuss the 
direct detection of dark matter. Since Ni are Majo- 
rana fermions, effective interaction terms N^Ni and 
Nicr^Ni vanish. Therefore, it's very hard to detect the 
elastic scattering between the target nuclei and dark 
matter Ni. One possibility is that the masses of N\ and 
N2 are highly degenerate. In this case, probably we can 
detect inelastic scatterings between Ni and N2. This 
case is discussed in Ref. [32]. Even in such a specific 
condition, the constraints given by direct detectors are 
still very limited. 

For further study, there are other potentially interest- 
ing phenomena including indirect signal from coannihi- 
lation and collider signature. For coannihilation events, 
data fitting is far more complicated due to the back- 
ground analysis. PAMELA collaboration reported their 
result indicating positron excess [33]. Ref. Q try to 
explain the excess within this model. But it needs an 
extremely large boost factor. For collider physics, some 
collider signatures of the model has been deduced in 
Ref. @,|2g AsLHC runs and collects data, this kind of 
events can be examined in the future. 
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